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High intensity inelastic neutron scattering experiments on the S—l quasi-one-dimensional bond- 
alternating antiferromagnet Ni(CgD24N4)(N02)C104 are performed in magnetic fields of up to 
14.8 T. Excitation in the high field magnetized quantum spin sohd (ordered) phase are investi- 
gated. In addition to the previously observed coherent long-lived gap excitation [M. Hagiwara et 
al., Phys. Rev. Lett 94, 177202 (2005)], a broad continuum is detected at lower energies. This 
observation is consistent with recent numerical studies, and helps explain the suppression of the 
lowest-energy gap mode in the magnetized state of NTENP. Yet another new feature of the exci- 
tation spectrum is found at slightly higher energies, and appears to be some kind of multi-magnon 
state. 



PACS numbers: 75.40.Gb, 75.10.Jm, 75.50.Ee 

I. INTRODUCTION 

The spin dynamics of gapped one-dimensional (1-D) 
quantum magnets in high magnetic fields has recently 
become an area of intensive research. The lowest-energy 
excitations in these systems are typically an S' = 1 triplet, 
and therefore subject to Zeeman splitting by an external 
field. With increasing field the gap in one of the three 
modes decreases and eventually vanishes, leading to a 1-D 
Bose condensation of magnonspi*Sii4 A qualitatively new 
state emerges at higher fields. In the presence of mag- 
netic anisotropy it is characterized by long-range Neel or- 
der. Despite that, its dynamical properties are nothing 
like those of conventional 3-D magnets. Quasiclassical 
spin wave theory (SWT) breaks down on the qualitative 
level, and the ground states are exotic "quantum spin 
solids" . Such is the case for the Haldane gap compounds 
NDMAF^i^i'^i^i^ and NENP,ia where aX H > there are 
three distinct magnon branches, while SWT can account 
for only two. 

The excitation spectrum in the high-field phase is 
non-universal, but can qualitatively vary from one spin 
system to another. In a recent work we investigated 
the bond-alternating dimerized S = 1 chain compound 
NTENPii-i2iia Like the structurally similar NDMAP, 
at = this material has a singlet ground state and 
a triplet of coherent gap excitations. It goes through 
a soft mode antiferromagnetic ordering transition at 
He ~ 9.7 T in a field applied parallel to the spin 
chains, and i/c ~ 13 T in a transverse fieldii However, 
its dynamical properties in a field are markedly differ- 



ent: only one coherent mode survives in the high-field 
phase of NTENP.^'^"* The highest-energy gap mode dis- 
appears well below He. This is due to interactions be- 
tween magnons that mix single- magnon and two-magnon 
states, opening an effective decay channel for the upper 
j-^Q^gji5ji6ji^ In addition, the lowest-energy mode (the 
one that goes soft), is a sharp long-lived excitation at all 
fields below the transition, but is absent in the magne- 
tized phase. Thus, only the middle mode persists above 
He. This one- mode picture is at present poorly under- 
stood. In the present paper we report a high intensity 
inelastic neutron scattering study that sheds new light 
on this peculiar behavior. At high fields, apart from the 
previously seen single coherent mode, we find a novel 
gapped low-energy excitation continuum. 



II. EXPERIMENTAL PROCEDURES AND 
RESULTS 

The structural elements of NTENP are discussed 
elsewhere, so only the key features are summarized 
here. NTENP crystallizes in the triclinic system 
(space group PI) with lattice constants a=10.747(l)A, 
&=9.413(2)A, c=8. 789(2)1, a=95.52(2)°, /3=108.98(3)° 
and 7=106.83(3)°. The ions are bridged by nitrito 
groups along the a axis with alternating bond distances 
of 2.142(3) and 2.432(6) A, and corresponding exchange 
constants are estimated as Ji = 2.1 meV and J2 = 
4.7 meV. Single-ion easy-plane anistropy of type DS^ 
is significant with D/J = 0.25, where J = (Ji -I- J2)/2 
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FIG. 1: A typical inelastic scan meaured at the ID AF zone- 
center in NTENP at ii" = 13.8 T> He (open circles). The 
background was measured away from the zone-center posi- 
tion (solid circles). Lines and shaded areas are guides for the 
eye. Inset: evolution of the FWHM resolution ellipsoid in the 
course of a typical inelastic scan used in this study. The solid 
line is the dispersion relation for the central coherent mode in 
NTENP at H = 13.8 T. The dashed hue represents the slope 
defined by the spin wave velocity. 



is the average coupling constant. In the experiments we 
employed a fully deuterated NTENP single crystal of ap- 
proximate volume 10 X 8 X 4 mm^. Sample environment 
was a 15 T split coil superconducting magnet with a dilu- 
tion refrigerator insert. Inelastic neutron data were col- 
lected at the IN14 cold-neutron spectrometer at Institut 
Laue-Langevin. In most cases neutrons with a fixed final 
energy (wave number) of 2.74 meV (1.15 A~^) were used 
with a Be higher-order filter positioned after the sam- 
ple and a horizontally focusing Pyrolitic Graphite (PG) 
analyzer. No additional beam coUimation devices were 
used. The sample was in all cases mounted with the b 
axis vertical (parallel to the applied field), allowing ac- 
cess to wave vectors in the (h, 0, 1) reciprocal-space plane. 

The data were collected in constant-^/ scans near the 
ID AF zone-centers qo = (1,0,0.2) and qi ^ (1,0,-2.2). 
Typical raw data are shown in Fig. ^ (open circles). 
The background (Fig. solid circles) was measured at 
q = (0.8,0,0.2) and q = (1.2,0,0.2), and found to be 
field-independent. A 5-point smoothed average of all 
background scans was subtracted from the data at each 
field. The resulting data sets collected at several fields 
are shown in Figs.|21and|31 They are to be compared with 
the data shown in Figs. 2 and 4 of Ref.lT3. that were mea- 
sured on a considerably less intense neutron source, with 
a somewhat broader energy resolution. 

The new data ioi H = 11.8 T ~ iJc are consistent with 
previous measurements. They clearly show a gap mode 
at fio; = 1.3 meV and a strong gapless feature at low en- 
ergies. The latter corresponds to the excitation branch 
that went soft at the critical field. At higher fields, the 
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FIG. 2: Background-subtracted inelastic scans measured at 
the ID AF zone-center qo = (1,0,0.2) in NTENP at various 
fields (open circles). Solid lines are fits to the data as de- 
scribed in the text. The shaded areas are the corresponding 
contributions of the three terms in the model cross section. 
The dashed line in (b)-(d) represent the experimental reso- 
lution, shown to emphasize the intrinsic broadening of the 
lower-energy peak. Arrows indicate the three main features 
of the spectrum. 
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FIG. 3: Background-subtracted inelastic scan measured at the 
ID AF zone-center qi = (1, 0, -2.2) in NTENP a± H = 14.8 T 
(open circles). Solid lines are fits to the data as described in 
the text. The shaded areas are the corresponding contribu- 
tions of the three terms in the model cross section. Arrows 
indicate the three main features of the spectrum. 
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scattering intensity at low energies decreases, but does 
not vanish entirely, as previously assumedii^ Instead, it 
is transformed into a broad peak (lower light shaded area 
in Fig. that is positioned below the strong coherent 
mode (dark shaded area in Fig. ^). The low intensity 
of this new feature explains why it could not be reliably 
detected in previous experiments: only a hint of a low- 
energy peak is visible in Fig. 4c of Ref. The new 
feature persists at higher fields, and its gap increases, 
with an intensity maximum shifting to around 0.6 meV 
at H = 14.8 T. At this field it is also observed at an 
equivalent but larger and almost orthogonal scattering 
vector qi — (1,0,-2.2) (Fig. Here the intensity is 
considerably smaller, consistent with its magnetic (as op- 
posed to vibrational) origin. It is important to emphasize 
that, unlike the sharp and intense lowest-energy mode in 
NDMAP, the low-energy feature in NTENP remains very 
broad and weak in the entire accessible field range. 

At the highest fields, yet another new feature of the 
spectrum is revealed aA, huj ~ 2 meV energy transfer (up- 
per light shaded area in Fig.^, and Figs.|2t,d). The new 
peak is also detected at qi = (1,0,-2.2) (Fig. EJ. This 
peak clearly does not correspond to the highest-energy 
member of the triplet seen at low fields. Indeed, the lat- 
ter vanishes well below He, and it energy extrapolates to 
much higher values with increasing field (somewhere in 
the range 3-3.5 meV at 14.8 T). Although at present we 
can not irrefutably prove that the new feature does not 
originate from spurious scattering, that possibility seems 
rather unlikely. Indeed, the peak is observed at several 
wave vectors, does not increase at a larger wave vector, 
shows a clear field dependence, and, as will be discussed 
below, appears to have a meaningful polarization char- 
acter. In the discussion below we shall assume that it is 
in fact real and magnetic in nature. 



III. DATA ANALYSIS AND DISCUSSION 

In order to learn more about the newly observed ex- 
citations, we performed fits to the experimental data, 
taking into account the effects of instrument resolution. 
The model cross section used for these fits included three 
separate branches, the partial contribution of each mode 
written as: 



da 



dfldE' 



1 



^2 , 



■L{1tIlI ~ hlUci,a,Ta) , (1) 

l+v^sin^qsi), (2) 



where {fiLOq^a) is the dispersion relation for mode a, 
L(x,T) is a normalized peak shape function of width 
r, taken in the Lorentzian or Dirac's ^-function (for 
r = 0) form, and v — 8.6 meV is the previously 
determined spin wave velocity4Siiii4 The model cross 
section was numerically convoluted with the resolution 
function of the spectrometer that was calculated in the 
Popovici approximation.^^ The intensities and widths of 



each mode, as well as the gap energies were refined to 
best-fit the experimental data at each field. The results of 
the fit are shown in heavy solid lines in Fig. 12 The grayed 
areas represent partial contributions of each mode. The 
excitation energies are plotted against the strength of ap- 
plied field in Fig. 01 (squares). Here we also show the data 
from previous studies (circles) lii^ 

The fits consistently yielded a zero intrinsic energy 
width for the central coherent mode. However, at H > 
He the lowest-energy feature is clearly broader than ex- 
perimental energy resolution. Its intrinsic energy full 
width at half height as about 0.30(6) meV at iJ = 12.8 T, 
and decreases to 0.15(4) meV at the highest accessible 
field of 14.8 T. The broad nature of this excitation is 
emphasized by the dashed lines in Fig. |21 that represent 
the calculated resolution width. We thus conclude that 
at H > He the lowest-energy excitations in NTENP are 
actually a continuum of states, rather than a single long- 
lived coherent mode as, for instance, in NDMAP. Such 
behavior is fully consistent with the Lanczos diagonal- 
ization results of Ref. Indeed, that numerical study 
predicted three coherent modes for NDMAP and only 
one for NTENP, but also a low-lying continuum of ex- 
citations in the latter system. In NTENP the observed 
intensity maximum of the low-energy spectral feature can 
be associated with the singularity on the lower bound of 
the continuum. This interpretation is consistent with the 
observed continuum spectrum being gapped or, at least, 
showing a pseudogap behavior. 

For the "2 meV" excitation, a zero intrinsic width 
was assumed for H = 13.8 T, where it is plainly too 
weak for a meaningful peak-shape analysis. However, at 
qi = (1,0,-2.2) at H = 14.8 T (Fig. EJ the feature is 
clearly broad beyond the effect of experimental resolu- 
tion. At this field its intrinsic width was determined to 
be 0.12(6) meV, comparable to that of the low-energy 
continuum. 

Additional information on the excitations in the high- 
filed phase was obtained in a partial polarization analysis 
at the highest accessible filed of iJ = 14.8 T. This pro- 
cedure involved a comparison of the inelastic intensities 
measured at qo = (1,0,0.2) (Fig. [21:) and the equivalent 
wave vector qi = (1,0,-2.2) (Fig. O. Note that, due 
to the large values of the crystallographic angle /3, these 
two wave vectors are almost orthogonal to each other, 
and form angle of 33° and 67° with the a (chain) axis, 
respectively. Neutron spectroscopy probes only the fluc- 
tuations of spin components perpendicular to the scat- 
tering vector. For qo the intensity scaling factors for ex- 
citation polarized along the x (parallel to the chains), y 
(in the horizontal plain, perpendicular to the chains) and 
z (vertical) axes are 0.30, 0.70 and 1.0, respectively. For 
qi the corresponding coefficients are 0.85, 0.15 and 1.0. 
Fits to the experimental data reveal that, corrected for 
resolution effects and the Ni^+ magnetic form factor, the 
intensities of all three components of the spectrum vary 
between the two wave vectors. At qi the corrected in- 
tensities of the central coherent mode and the low-energy 
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FIG. 4: Field dependence of excitation energies at the ID 
AF zone-center in NTENP. Open squares are positions of 
the intensity maximum in the low-lying excitation continuum 
(grayed area). Gray squares are the energies of the coherent 
central mode. Solid squares correspond to the "2 meV" peak. 
Circles are from Ref. ll4L 

continuum are decreased by factors of 1.6(1) and 1.27(4), 
respectively, as compared to qo. In contrast, after form 
factor and focusing corrections, the "2meV" excitation 
is stronger at qi, by a factor of 2.4(3). This intensity 
pattern is consistent with the "2meV" excitation being 
polarized parallel to the spin chains, and the two lower- 
energy modes being transverse in nature. 

Regarding the origin of the "2 meV" feature, we note 



that its peak energy is approximately equal to the sum 
of the gap in the main coherent mode and the energy of 
the intensity maximum of the continuum. This may indi- 
cate that the "2 meV" mode is a composite excitation of 
these two lower-energy states, similar to bound magnon 
states in the classical magnet TMMC.^2, Such excitations 
near q — n are permitted by energy-momentum conserva- 
tion laws. Indeed, the intrinsic breaking of translational 
symmetry (bond-alternation) in NTENP makes the wave 
vectors q = and q = n equivalent, allowing excitations 
of total momentum tt composed of a continuum state 
near q — and a single-magnon state at g = tt, or vice 
versa. A further theoretical study is required to clarify 
this point. 

IV. CONCLUSION 

The examples of NDMAP and NENP may lead to a 
simple-minded conclusion that the high field phase of 
anisotropic gapped quantum spin chains is entirely dom- 
inated by long-lived coherent excitations. The new re- 
sults on NTENP clearly demonstrate the opposite: broad 
spectral continua may, in fact, play a key role in the dy- 
namics of highly anisotropic quantum spin solids even at 
low energies. 



Acknowledgments 

We thank Sei-Ichiro Suga and Takahumi Suzuki for 
fruitful discussions. This work was in part supported by 
the Molecular Ensemble research program from RIKEN 
and the Grant-in- Aid for Scientific Research on Priority 
Areas from the Japanese Ministry of Education, Culture, 
Sports, Science and Technology. Work at ORNL was 
sponsored by the Division of Materials Sciences and En- 
gineering, Office of Basic Energy Sciences, U.S. Depart- 
ment of Energy, under contract DE-AC05-00OR22725 
with Oak Ridge National Laboratory, managed and op- 
erated by UT-Battelle, LLC. 



1 H. J. Schulz, Phys. Rev. B 34, 6372 (1986). 

^ K. Katsumata, H. Hori, M. D. T. Takeuchi, A. Yamagishi, 

and J. R Renard, Phys. Rev. Lett. 63, 86 (1989). 
^ I. Affleck, Phys. Rev. Lett. 65, 2477 (1990). 

* M. Takahashi and T. Sakai, J. Phys. Soc. Jpn. 60, 760 
(1991). 

^ Z. Honda, H. Asakawa, and K. Katsumata, Phys. Rev. 
Lett. 81, 2566 (1998). 

® Y. Chen, Z. Honda, A. Zheludev, C. Broholm, K. Kat- 
sumata, and S. M. Shapiro, Phys. Rev. Lett. 86, 1618 
(2001). 

A. Zheludev, Y. C. Z. Honda, C. Broholm, and K. Kat- 
sumata, Phys. Rev. Lett. 88, 077206 (2002). 

* A. Zheludev, Z. Honda, C. Broholm, K. Katsumata, S. M. 
Shapiro, A. Kolezhuk, S. Park, and Y. Qiu, Phys. Rev. B 



68, 134438 (2003). 

^ A. Zheludev, cond-mat/0507534 (2005). 

" L.P. Regnault, I. Zaliznyak, J. P. Renard, and C. Vettier, 

Phys. Rev. B 50, 9174 (1994). 
^ Y. Narumi, M. Hagiwara, M. Kohno, and K. Kindo, Phys. 

Rev. Lett. 86, 324 (2001). 
^ A. Zheludev, T. Masuda, B. Sales, D. Mandrus, T. Papen- 

brock, T. Barnes, and S. Park, Phys. Rev. B 69, 144417 

(2004). 

^ L. P. Regnauh, M. Hagiwara, N. Metoki, Y. Koike, K. 
Kakurai, and A. Stunauh, Physica B 66, 350B (2004). 

^ M. Hagiwara, L. P. Regnault, A. Zheludev, A. Stunault, 
N. Metoki, T. Suzuki, S. Suga, K. Kakurai, Y. Koike, 
P. Vorderwisch, et al., Phys. Rev. Lett. 94, 177202 (2005). 

^ T. Suzuki and S. Suga, Phys. Rev. B 72, 014434 (2005). 



5 



A. Kolozhuk and S. Sachdov, cond-mat/0511353. 
M. E. Zhitomirsky, cond-mat/0601405. 
M. Popovici, Acta Cryst. A31, 507 (1975). 



I. Heilmann, J. Kjoms, Y. Endoh, G. Reiter, G. Shirane, 
and R. Birgeneau, Phys. Rev. B 24, 3939 (1981). 



